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I.  INTRODUCTION 


Several  previous  Desmatics  technical  reports  have  Investigated  the 
use  of  a  logistic  function  in  the  development  of  indirect  impact  acceler¬ 
ation  injury  prediction  models.  Empirical  data  have  been  used  to  fit  these 
models  and  evaluate  their  ability  to  accurately  predict  fatal  injuries.  The 
most  recent  of  these  reports  [3]  explored  the  possibility  of  developing  a 
model  based  on  head  dynamic  response  data.  Models  using  only  sled  acceler¬ 
ation  profile  terms  were  also  considered ,  as  were  models  using  a  mixture  of 
the  two  types  of  variables.  The  results  of  that  study  have  been  influential 
in  the  selection  of  variables  to  be  considered  in  this  investigation. 

The  prediction  models  used  in  this  report  are  of  the  same  functional 
type  as  considered  previously: 

k  -1 

P(x)  «  {  1  +  exp [  -  (  3n  +  £  )]  } 

u  i*l  x  1 

where : 

x  -  (x^.x^, . . . .x^)  denotes  the  set  of  Independent  variables, 

• • »8r)  denotes  a  set  of  parameter  values, 
and  P(x)  denotes  the  true  probability  of  Injury  corresponding  to  x. 

A  discussion  of  how  this  type  of  model  may  be  related  to  the  forces  and 
moments  in  the  region  of  injury  may  be  found  in  [1]. 

The  generic  type  of  model  described  above  has  previously  been  applied 
to  a  set  of  63  -C  accelerator  runs  involving  Rhesus  monkeys  with  securely 
restrained  torsos  and  unrestrained  heads.  The  data  were  collected  by  the 
Naval  Biodynamics  Laboratory  (NBDL)  as  part  of  its  research  effort  on  im¬ 
pact  acceleration  injury  prevention.  Three  sled  profile  variables  and  twelve 

head  dynamic  response  variables  were  considered  for  inclusion  in  the  models. 
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In  the  models  which  included  both  types  of  variables,  peak  sled  acceleration 
was  found  to  be  the  most  Important  variable  while  the  peak  z- component  of 
head  angular  velocity  also  contributed  to  the  explanatory  power  of  the  model. 
When  only  head  dynamic  response  variables  were  used,  the  peak  y -component  of 
head  angular  acceleration  was  also  found  to  be  important. 

Research  to  date  has  found  that  peak  sled  acceleration  is  a  better  pre¬ 
dictor  of  fatality  in  -Gx  Rhesus  runs  than  is  any  combination  of  head  dynamic 
response  variables  yet  considered.  Furthermore,  inertial  data  are  not  available 
for  several  high  level  (>85  G)  experiments.  Since  the  most  useful  data  for  the 
development  of  a  predictive  model  are  those  obtained  near  the  boundary  between 
injury  and  noninjury,  the  high  level  experiments  are  of  primary  interest  in 
the  formulation  of  the  logistic  models  considered  here.  This  combination  of 
factors  leads  to  the  consideration  of  models  which  do  not  rely  on  Inertial 
data. 

Although  head  dynamic  response  data  are  not  available  for  many  of  the 
accelerator  runs  of  greatest  interest,  there  are  photographic  data  available 
from  which  the  initial  yaw  and  pitch  angles  of  the  head  may  be  determined. 

As  mentioned  in  [1] ,  the  initial  yaw  angle  should  be  highly  correlated  with 
the  z-component  of  head  angular  velocity.  It  also  seems  reasonable  that 
the  initial  pitch  angle  might  be  correlated  with  the  y-component  of  angular 
acceleration.  Thus,  inclusion  of  these  initial  head  conditions  should  at 
least  partially  compensate  for  the  lack  of  inertial  data. 

The  data  base  used  in  this  Investigation  consists  of  23  high  level  -G 

experiments.  For  Run  No.  LX0657  the  initial  pitch  angle  is  unavailable 

and  for  Run  No.  LX3010  the  initial  yaw  angle  is  not  known.  Thus,  only  21 

data  vectors  can  be  used  in  the  analysis.  Thirteen  of  these  sled  runs  have 

been  used  in  previous  analyses  while  the  remaining  eight  experiments  could 
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not  be  used  previously  owing  to  the  lack  of  inertial  data.  The  independent 
variables  used  in  this  analysis  are  peak  sled  acceleration  (G)»  initial  yaw 
angle  of  the  head  (deg),  and  initial  pitch  angle  of  the  head  (deg).  Because 
of  symmetry  about  the  mid-sagittal  plane  in  -G  accelerator  runs,  the  absolute 
value  of  the  initial  yaw  angle  was  the  variable  actually  used  in  the  models. 

As  in  previous  studies,  the  dependent  variable  is  binary  (1 »  fatal  run,  0* 
non-fatal  run) .  A  listing  of  the  data  set  is  given  in  Table  1. 
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Run  Subject  Observed  ^  ^  Peak  Sled  Initial  Initial 

Number  Number  Probability  P(x, ,x_)  P(x, ,x„,x„)  Acceleration  Yaw  Pitch 
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Table  1:  The  Data  Set  With  Predicted  Probabilities  of  Fatality  For  Two  Models 


II.  MODEL  CONSTRUCTION 


As  mentioned  in  the  previous  section ,  the  data  base  used  in  this  in¬ 
vestigation  consists  of  21  complete  observations.  Only  16  animals  mere 
used  for  this  set  of  experiments,  so  some  dependence  exists  in  the  data  as 
a  result  of  having  multiple  sled  runs  with  the  same  animal.  As  stated  in 
[3],  this  type  of  dependence  will  tend  to  make  the  model  conservative;  i.e., 
the  probability  of  injury  will  be  overpredicted  for  a  specified  x. 

Initially,  a  model  was  considered  which  included  the  three  independent 
variables  and  all  of  their  two-way  interactions.  A  backward  selection  method 
was  employed  in  order  to  eliminate  those  terms  which  did  not  contribute  sig¬ 
nificantly  to  the  model's  ability  to  explain  the  occurrence  or  nonoccurrence 
of  fatal  injury.  In  this  way,  all  of  the  interactions  were  found  to  be  un¬ 
important  terms  in  the  model. 

Denote  the  independent  variables  by  x^,  x^,  and  x^,  where: 

x^«peak  sled  acceleration, 

x^  ■  absolute  value  of  the  initial  yaw  angle  of  the  head, 
and  x^ *  initial  pitch  angle  of  the  head. 

The  predicted  probability  of  a  fatality  based  on  these  three  variables  is: 

£(x^,  X£,  x^)  *  {l  +  exp[-(-18.4  +  .157  x^  +  .0352  -  .0222  x^)  ]  }  ^  • 

A  comparison  between  observed  probabilities  and  predicted  probabilities  for 
this  model  is  given  in  Table  1. 

The  variable  which  has  the  least  influence  on  ?(x^,  x^,  x^)  is  x^  ,  the 
Initial  pitch  angle.  A  test  of  whether  x^  is  making  a  significant  contribution 
to  the  model  can  be  performed  by  calculating  the  following  quantities: 


-  -2  log  likelihood  for  model  containing  (x^,  x2) 
and  -  -2  log  likelihood  for  model  containing  (Xj,  x2»  x^). 

Under  the  null  hypothesis  that  x3  does  not  improve  the  model,  the  statistic 
-L2  has  an  approximate  Chi-square  distribution  vith  1  degree  of  freedom. 

For  this  data  set,  the  statistics  are  given  below: 

1^-17.695  ,  L2«  17.213,  1^-12-0.482  . 

Therefore,  there  is  no  significant  effect  due  to  x3  (p  -  .487),  given  that  xx 
and  x2  are  already  in  the  model. 

Since  x3  does  not  significantly  improve  the  model's  ability  to  explain 
the  data,  it  is  probably  better  to  use  a  model  based  only  on  x^  and  x2  .  The 
predicted  probability  of  a  fatality  for  this  reduced  model  is: 

f(xlt  x2)  -  {  1  +  exp [-(-18.5+  .152  x^.0442  x^]}"1. 

The  predicted  probabilities  for  this  model  are  also  given  in  Table  1.  Both 
and  x2  are  important  in  this  model  so  no  further  reduction  is  possible. 

The  two  models  given  in  this  section  explain  the  variation  in  the  data 
about  equally  well.  However,  when  using  these  models  to  predict  new  obser¬ 
vations,  the  simpler  model  is  probably  more  trustworthy.  The  small  effect 
attributable  to  the  pitch  angle  of  the  head,  given  that  peak  sled  accelera¬ 
tion  and  initial  yaw  angle  have  already  been  taken  into  account,  could  easily 
be  spurious.  While  any  variable  added  to  a  model  will  increase  the  ability 
of  that  model  to  explain  the  data,  the  new  variable  might  be  useless  as  a 
predictor  of  future  observations.  When  the  perceived  effect  of  a  factor  is 
very  small,  that  effect  is  more  likely  to  be  the  result  of  random  confound¬ 
ing  of  that  factor  with  some  variable  which  has  not  been  Included  in  the  anal¬ 
ysis.  The  confounding  can  result  in  an  extremely  misleading  prediction  model. 
Therefore,  of  the  two  models  given  here,  the  two-variable  model  is  recoonended. 
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Attention  will  be  restricted  to  it  in  the  remainder  o£  this  report. 

The  predicted  probability  of  a  fatality  under  the  specified  model 
is  as  given  above,  where  x^  is  the  peak  sled  acceleration  and 

x ^  ie  the  absolute  value  of  the  initial  yaw  angle  of  the  head.  The  pre¬ 
dicted  region  of  fatality  is  the  half  plane  defined  by: 

.152  x^  +  .0442  x^  c  . 

The  constant  c  is  determined  by  the  particular  predicted  probability  one 
wishes  to  use  as  a  cutoff  value.  For  example,  *  .5  corresponds 

to  c  *  18.5  . 
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III.  COMPARISON  WITH  EARLIER  MODELS 


In  this  section  a  comparison  Is  made  between  the  model  developed 

in  this  report  and  two  of  the  models  presented  in  [3].  As  an  aid  in  this 

comparison,  it  has  been  found  useful  to  estimate  one  other  model.  This 

last  model  includes  only  peak  sled  acceleration  as  an  independent  variable 

but  is  based  on  the  entire  set  of  available  -G  accelerator  data.  As  be- 

x 

fore,  a  logistic  response  model  is  assumed  to  be  appropriate. 

Since  the  models  developed  in  [3]  were  based  on  a  different  set  of 
data,  they  cannot  be  compared  directly  with  the  model  presented  in  Section 
II.  However,  there  are  thirteen  observations  common  to  the  two  data  bases 
and  some  assessment  of  relative  performance  may  be  obtained  by  comparing  the 
predicted  probabilities,  using  the  various  models,  for  this  subset  of  the 
data.  In  order  to  facilitate  this  comparison,  the  following  notation  is 
defined: 

y^  ■  absolute  value  of  the  peak  z-component  of  head  linear  acceleration, 

■  peak  head  linear  resultant  acceleration, 
y^ *  absolute  value  of  the  peak  y-component  of  head  angular  acceleration, 
and  y^ » absolute  value  of  the  peak  z-component  of  head  angular  velocity. 

As  before,  and  x^  denote  peak  sled  acceleration  and  the  absolute  value 
of  the  initial  yaw  angle  of  the  head,  respectively. 

In  keeping  with  the  notation  used  earlier  in  this  report,  let  P(x^,x2) 
denote  the  predicted  probability  of  a  fatality  for  the  model  given  in  Section 
II.  Further,  consider  the  best  three-variable  head  dynamic  response  model, 
given  in  [3],  which  involves  the  three  variables  y^,  y 2»  and  y^.  Denote 

A 

the  predicted  probability  of  a  fatality  using  this  model  by  Piyj.yj.y-j). 


Also  consider  the  best  two-variable  mixed  model,  given  in  [3],  which  in¬ 
volved  and  y^.  Let  P(x^,y^)  denote  the  predicted  probability  for  this 
model.  A  comparison  of  these  predicted  probabilities  for  the  thirteen 
common  data  vectors  is  given  in  Table  2.  (Minor  inaccuracies  have  been 
found  in  the  earlier  technical  report.  The  predicted  probabilities  given 
here  are  therefore  slightly  different  from  those  given  previously.) 

None  of  the  three  models  compared  in  Table  2  is  a  better  predictor 
of  fatality  for  each  of  the  data  points  under  consideration.  It  is  there¬ 
fore  necessary  to  use  some  form  of  summary  statistic  to  aid  in  the  compari¬ 
son.  Perhaps  the  most  obvious  measure  of  performance  is  given  by: 

13  ~  2 
Q.  -  E<0  -p  r , 

1  i-1  1  1 

where  0i  is  the  observed  probability  of  a  fatality  (0  or  1)  for  observation 
i  and  P^  is  the  corresponding  predicted  probability.  However,  use  of 
has  one  serious  drawback.  Observations  with  predicted  probabilities  near 
.5  automatically  contribute  large  terms  to  the  sum,  since  the  observed  prob¬ 
ability  must  be  either  0  or  1.  It  seems  reasonable  that  the  statistic  used 
should  only  be  large  if  an  observed  fatal  run  has  a  predicted  probability 
less  than  .5  or  an  observed  nonfatal  run  has  a  predicted  probability  greater 
than  .5  .  One  statistic  of  this  type  is: 

13  (0  -P  )2 

Q  -  I  -i - i -  * 

i-lPi(l-?i) 

The  weights  used  in  this  sum  are  the  inverse  variances  of  binomial  observa- 

A 

tions  where  is  the  probability  that  the  observation  is  1. 

Q.j  and  have  been  calculated  for  the  three  models  being  compared.  The 
corresponding  values  are: 
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Run 

Number 

Subject 

Number 

Observed 

Probability 

^(y1.y2*y3) 

^<x1,y2) 

LX3188 

AR8863 

0 

.1179 

.3241 

.2022 

LX1362 

A03935 

0 

.1190 

.3889 

.3175 

LX1894 

A03933 

0 

.6328 

.4996 

.3175 

LX1359 

A04099 

0 

.1620 

.3326 

.3866 

LX1893 

A03924 

0 

.3289 

.5324 

.3866 

LX3709 

AR3790 

1 

.2116 

.0800 

.0776 

LX3192 

AR8866 

1 

.2562 

.3180 

.5265 

LX1363 

A03935 

1 

.0845 

.7345 

.5463 

LX1905 

A04101 

1 

.9450 

.8081 

.6550 

LX1896 

A03946 

1 

.6962 

.9028 

.8023 

LX1365 

AO 3921 

1 

1.0000 

.8699 

.8850 

LX1895 

A03951 

1 

.9999 

.9637 

.9954 

LX1360 

AO4099 

1 

.9106 

.9628 

.9928 

Table  2:  Comparison  of  Predicted  Probabilities  for  Observations 
Common  to  the  Two  Data  Sets. 


Model 


2 


yl*y2,y3  2,679  20,732 

xx,y4  2.348  18.328 

xltx2  1.994  16.976 


Both  and  measure ,  In  some  sense,  lack  of  fit  of  the  model  being  con¬ 
sidered.  That  is,  they  are  large  when  the  model  is  not  adequately  explain¬ 
ing  the  observed  pattern  of  fatalities.  Thus,  using  either  criterion. 


P(x^,x^)  performs  better  for  this  subset  of  the  data  than  does  either  of 
the  models  taken  from  [3].  Also,  Ptx^.y^)  explains  the  data  better  than 
does  ?(y1,y2,y3) ,  as  could  be  expected  from  previous  results. 

A 

There  are  at  least  two  possible  reasons  for  the  fact  that  P(x^,x2) 
performs  better  than  P(x^,y4).  First,  it  is  possible  that  the  initial  yaw 
angle  of  the  head  is  a  better  predictor  of  injury  than  is  the  peak  z-com- 


ponent  of  head  angular  velocity.  However,  it  is  also  possible  that  the 
peak  z-component  of  head  angular  velocity  is  the  better  predictor  but  that 
its  coefficient  in  the  model  is  not  as  well  estimated.  This  circumstance 
could  easily  arise  as  a  consequence  of  the  fact  that  different  data  sets 
were  used  to  estimate  the  two  models. 

A 

The  data  set  used  to  estimate  the  coefficients  in  P(x^,y4>  consisted 

A 

of  63  observations  while  only  21  observations  were  used  to  estimate  P(x^,x2>. 

However,  those  21  observations  were  all  high-level  experiments  and  thus  nearer 

to  the  boundary  between  injury  and  non-injury.  That  type  of  observation  is 

more  useful  for  estimating  the  coefficients  in  a  logistic  response  model. 

It  is  therefore  necessary  to  quantify  the  usefulness  of  each  observation  in 

order  to  determine  the  relative  value  of  the  two  data  sets.  This  is  usually 

accomplished  by  calculating  some  function  of  Fisher's  information  matrix  [2]. 
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Suppose  w  is  the  dependent  variable  to  be  measured.  Then  for  any 
particular  set  of  values  of  the  independent  variables  x,  w  has  a  probability 
distribution  which  may  depend  on  the  vector  of  parameters  ji.  Denote  the 
probability  density  function  of  w  at  x  by  p(w|x»£)  .  For  a  single  observa¬ 
tion  on  w  at  x,  the  Fisher  information  matrix  for  j}  is  that  matrix  whose 
(i+l,j+l)  element  is  given  by: 

E  {-  32  log  p(w|x,J)/36136j}  , 

where  ••• »Bk)  and  E  denotes  the  expectation  over  the  probability 

distribution  of  w.  The  information  matrix  of  N  Independent  observations 

N 

taken  respectively  at  * £2  *  * '  *  ■  IJ(x^;jJ),  where  J ;J3)  is 

the  information  matrix  for  a  single  observation,  as  given  above,  and 
X  * (x^,x^,» .. ,x^)  .  A  more  complete  discussion  of  this  concept  may  be  found 
in  [2]. 

For  the  injury  prediction  models  being  considered  here,  w  is  binary 
(1  ■  fatality,  0  -  nonfatality) .  The  probability  density  function  of  w, 
given  x,  is: 

p(w|x,B)  -  [P(x)]W[l-P(x)]1~W;  w-0,1  , 
k  -1 

where  P(x)  -  {  l  +  expt~Bn+  IB.x  )]  } 

i«i  1 

For  this  model,  with  Xq«1,  the  (i+l,j+l)  element  of  J(x;S)  is: 

xixjp(*Hl  -  P(x)  1  . 

Since  J(x;ji)  depends  on  POO,  the  information  matrix  cannot  be  computed  unless 
the  correct  model  is  known.  It  is  possible,  however,  to  obtain  some  indica¬ 
tion  of  how  much  information  is  available  in  the  samples  by  using  the  estimate 
of  P(jc).  This  estimate  must,  of  course,  ba  applicable  to  both  data  sets  so 
that  only  the  estimated  information  matrices  will  be  directly  comparable. 

The  only  variable  which  is  available  for  all  of  the  observations  in  both 
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data  sets  is  peak  sled  acceleration.  Therefore,  the  estimate  of  P(x)  to 
be  used  in  computing  the  information  matrices  must  be  a  function  of  this 
variable  alone.  Fortunately,  peak  sled  acceleration  is  known  to  be  a  good 
predictor  of  fatality.  Furthermore,  93  total  observations  are  available 
for  the  estimation  procedure,  so  that  the  coefficients  in  the  model  should 
be  fairly  well  estimated. 

The  estimated  model  using  only  peak  sled  acceleration  as  an  independent 
variable  is: 

P(x)  -  {  l  +  exp(13.6-  .120x)}_1  . 

The  predicted  values  from  this  model  can  now  be  used  to  estimate  J(x;j3)  for 
each  of  the  93  observations.  Three  different  information  matrices  have  been 
calculated  and  are  given  below.  1^  refers  to  the  complete  set  of  93  obser¬ 
vations,  refers  to  the  62  observations  used  to  fit  models  involving  head 
dynamic  response  variables,  and  1^  refers  to  the  set  of  21  observations  used 
to  estimate  the  models  given  earlier  in  this  report. 


3.72 

411. 25\ 

411.25 

46458.52 / 

2.36 

259. 22\ 

259.22 

20019.05/ 

1  3.13 

356. 87\ 

,356.87 

41174.91 / 

The  inverse  of  the  information  matrix  provides  a  lower  bound  for  the 

variance  matrix  of  an  unbiased  estimate  of  j5.  It  is  therefore  desirable 

for  the  inverse  to  be  "small”  or,  equivalently,  for  the  information  matrix 

to  be  "large"  in  some  sense.  One  popular  optimality  criterion  is  D-opti- 

raality,  which  says  to  maximize  the  determinant  of  the  information  matrix. 

Since  |l^|  *1521  and  jl^l  *  1290,  the  sample  used  in  this  report  contains 
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more  information,  according  to  this  criterion,  than  does  the  sample  used 
in  the  previous  technical  report.  It  is  therefore  not  surprising  that  the 
model  developed  here  performs  better  on  the  subset  common  to  both  samples. 
However,  D-optimality  is  by  no  means  the  only  criterion  used  to  evaluate 
designs.  A  number  of  different  functions  of  the  information  matrix  can  be 
used.  The  only  time  one  matrix  can  be  said  to  be  larger  than  another  is 
when  their  difference  is  positive  definite.  In  this  case  l^-I^  is  not 
positive  definite  (the  determinant  is  negative)  so  it  is  impossible  to  say 
definitively  which  sample  contains  more  information. 


IV.  SUMMARY 


Using  a  data  base  of  21  high-level  accelerator  runs,  models  have  been 
considered  which  included  peak  sled  acceleration,  initial  yaw  angle  of  the 
head,  and  initial  pitch  angle  of  the  head  as  Independent  variables.  While 
initial  yaw  was  found  to  be  a  useful  predictor  of  Injury,  the  initial  pitch 
angle  was  found  to  make  little,  if  any,  additional  contribution  to  the  model, 
given  that  the  other  variables  had  already  been  included.  Furthermore,  no 
significant  interactions  between  the  three  variables  were  found  to  exist. 

In  a  previous  technical  report,  [3],  models  were  developed  which  de¬ 
pended  on  the  use  of  head  dynamic  response  data.  Of  the  63  data  vectors 
used  in  that  investigation,  thirteen  observations  Included  head  Initial 
conditions  as  well.  The  model  developed  in  this  report  was  compared  to  the 
earlier  models  on  the  basis  of  each  model's  ability  to  predict  injury  for 
the  thirteen  common  observations.  Using  two  different  relative  performance 
criteria,  the  initial  conditions  model  was  found  to  perform  better  than  the 
others. 

Since  the  model  developed  in  this  report  used  a  different  set  of  data 
than  the  earlier  models,  there  was  some  concern  that  its  superior  performance 
on  the  common  data  set  might  be  the  result  of  the  model  being  better  estimated, 
rather  than  reflect  the  relative  worth  of  the  variables.  Therefore,  an  at¬ 
tempt  was  made  to  quantify  the  amount  of  information  available  in  each  sample. 
Fisher's  information  matrix  was  calculated  for  each  sample,  using  an  estimate 
of  the  predicted  probability  for  each  data  point.  These  estimates  were  taken 
from  a  single-variable  model  using  peak  sled  acceleration  and  based  on  all 
available  -Gx  acceleration  data.  Unfortunately,  no  clear  decision  could  be 
made  as  to  which  sample  contains  the  most  information.  Using  at  least  one 


criterion,  however,  the  data  base  used  in  this  report  is  better.  Thus,  the 
model  that  performs  best  on  the  comnon  observations  is  the  model  based  on: 

-  peak  sled  acceleration 

and  -  absolute  value  of  initial  yaw  angle  of  the  head. 

However,  because  at  least  part  of  its  performance  may  be  due  to  more  Inherent 
information  content  in  the  observations  on  which  it  was  based,  there  is  no 
conclusive  evidence  that  x^  and  x^  are  the  best  predictors. 
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